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ABSTRACT 

We investigate the molecular environment of the Galactic supernova remnant (SNR) Kesteven 78 
and perform an XMM-Newton X-ray spectroscopic study for the northeastern edge of the remnant. 
SNR Kes 78 is found to interact with the molecular clouds (MCs) at a systemic local standard of rest 
velocity of 81 kms~^. At around this velocity, the SNR appears to contact a long molecular strip in the 
northeast and a large cloud in the east as revealed in the ^^CO line, which may be responsible for the 
radio brightness peak and the OH maser, respectively. The ^^CO-line bright region morphologically 
matches the eastern bright radio shell in general, and the SNR is consistent in extent with a CO 
cavity. Broadened "'^^CO line profiles discerned in the eastern maser region and the western clumpy 
molecular arc and the elevated ^^CO J=2-l/ J=l-0 ratios along the SNR boundary may be signatures 
of shock perturbation in the molecular gas. The SNR-MC association places the SNR at a kinematic 
distance of 4.8 kpc. The X-rays arising from the northeastern radio shell are emitted by underionized 
hot (~ 1.5 keV), low-density 0.1 cm~'^) plasma with solar abundance, and the plasma may be of 
intercloud origin. The age of the remnant is inferred to be about 6 kyr. The size of the molecular 
cavity in Kes 78 implies an initial mass around 22Mq for the progenitor. 

Subject headings: ISM: individual objects (G32.8— 0.1 — Kes 78), — ISM: molecules, — supernova 
remnants 



1. INTRODUCTION 

Born in giant molecular clouds (MCs), massive stars 
evolve from collapsing molecular cores to core-collapse 
supernovae. During their short lifetimes, they evacuate 
a cavity with their energetic winds and ionizing radia- 
tion, and later the supernova shocks will often collide 
with the cavity wall. In the past three decades, dozens 
of supernova remnants (SNRs) were found to be in phys- 
ical contact with MCs (see Jiang et al. 2010 and refer- 
ences therein). In particular, the 1720 MHz OH masers, 
collisionally excited by the shock-heated molecular gas, 
have been regarded as the signposts of SNR-MC interac- 
tion (e.g.. Frail et al. 1996; Lockett et al. 1999; Frail & 
Mitchell 1998; Wardle & Yusef-Zadeh 2002). Recently, 
several GeV and TeV sources were found to spatially cor- 
respond to MCs that are associated with SNRs, such as 
IC443 and W28, in which p — p collision plays an im- 
portant role in the gamma-ray emission. Detailed inves- 
tigation of the molecular environment of the interacting 
SNRs is therefore of great interests. 

Kesteven 78 (G32. 8-0.1), with a partly brightened ra- 
dio shell, was identified as a Galactic SNR by Caswell 
et al. (1975) based on the 408 and 5000MHz observa- 
tions. A distance of ~ 8 pc was estimated using the 
radio surface brightness - diameter (S-D) relationship 
(Caswell et al. 1983). A 21cm HI observation with a res- 
olution 2' X 130' X 6.3 km s^^ was made by Gosachinskii 
& Khersonskii (1985), which suggested a corresponding 
local standard of rest (LSR) velocity of 90kms~^, a dis- 
tance of 9 kpc, and an age of 12 kyr for the SNR. Kas- 
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sim (1992) derived the spectral index as —0.5. A single 
1720 MHz hydroxyl (OH) maser spot at LSR velocity 
86.1 kms~^ was detected on the radio shell by Koralesky 
et al. (1998) in the Very Large Array (VLA) observations, 
which indicates a shock interaction with MCs; a distance 
of 5.5 or 8.8 kpc was then suggested on the assumption 
that the shock propagation is perpendicular to line of 
sight (LOS). A classical antisystematic S-shaped profile 
of Stokes V was also detected, suggesting the presence of 
Zeeman splitting. A millimeter-wavelength ^^CO ( J=l- 
0) observation was carried out by Zhou et al. (2007), 
which suggests that an eastern molecular gas detected in 
the velocity range 72-88 kms~^ is related to the SNR, 
and no line broadening was found from it. In a Spitzer- 
IRAC mid-infrared (IR) survey of Galactic SNRs, which 
was made to search for IR counterparts with existing ra- 
dio images, Kes 78 is classified as "not detected but con- 
fused" one (Reach et al. 2006). In the optical band, fila- 
mentary and diffuse emission is detected in the SNR; the 
[SII]/Ha ratio higher than 1.2 is extensively found, sug- 
gesting that the optical emission arises from the shock- 
heated gas (Boumis et al. 2009, B09). An extended very 
high energy (VHE) 7-ray source HESS J1852-000 is re- 
vealed by the H.E.S.S team to be close to the east edge 
of the remnant 

Motivated by the indication of the OH maser for SNR- 
MC interaction; a hint of the previous CO observation 
for a related MC in the east; and the irregular morphol- 
ogy of the SNR, which implies a complex, non-uniform 
ambient medium, we examine the detailed distribution 
of the environmental molecular gas and its physical rela- 
tion with Kes 78, based on both our new and the archival 
observations of three CO rotational transitions and the 

■'http: / /www. mpi-hd.mpg.de/hfm/HESS/pages/home/som/2011/02/ 
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archival data of the XMM-Newton X-ray observation. In 
Section 2, we briefly describe the observations and data 
reduction. Our results are presented in Section 3. Physi- 
cal properties of the SNR are discussed in Section 4. The 
conclusions are summarized in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. CO Observations and Data 

The observations of millimeter molecular emissions 
toward SNR Kes 78 were first made in 2009 Novem- 
ber with the 13.7 m millimeter-wavelength telescope of 
the Purple Mountain Observatory at Delingha (here- 
after PMOD), China. A superconductor- insulator- 
superconductor (SIS) receiver and two acousto-optical 
spectrometers (AOSs) were used to simultaneously ob- 
serve the i^CO (J^l-^O) line (at 115.271 GHz) and the 
13C0 (J=l-0) line (at 110.201 GHz). We mapped the 
radio bright shell region of Kes 78 with 1' grid spacing 
and a 39' x 27' large region centered at (18''51'"03'.97, 
-00°12'22".5, J2000.0) with 2' grid spacing. The 
AOS bandwidth and velocity resolution are 145 MHz 
and 0.37kms-i for ^^CO{J^1-0) and 43 MHz and 
0.11 kms~^ for ^'^CO (J=l-0). During the observation 
epoch, the half-power beam width (HPBW) of the an- 
tenna was 56", the main beam efficiency was about 62% 
at zenith, and the pointing accuracy was better than 5". 
The typical system temperature was around 230 K. The 
observed LSR velocity ranges were —102 to 276kms~^ 
for i^CO (J=l-0) and 28 to 143 kms'^ for ^^CO (J=l- 
0). 

The follow-up observation was made in the ^^CO 
( J=2-l) line (at 230.538 GHz) during 2010 January and 
February using the Kolner Observatory for Submillime- 
ter Astronomy (KOSMA) 3m submillimeter telescope in 
Switzerland. An SIS receiver and a medium-resolution 
AOS spectrometer were used. A 31' x 31' area centered 
at (18^51™03".97, -00°08'22".5, J2000.0) was mapped 
with a grid spacing of 1' in on-the-fly mode. The HPBW 
of the telescope was 130", the main beam efficiency was 
54% during our observation, and the pointing accuracy 
was about 10". The AOS bandwidth and velocity reso- 
lution were about 300 MHz and 0.2kms~^. 

All the CO data were reduced with GILD AS/ CLASS 
package El and analyzed with IDL and KARMA (Gooch, 
1996). After the baseline subtraction and the calibration 
for main beam efficiency and elevation, the spectra were 
resampled to a uniform velocity resolution of 0.5kms~^ 
and the data were convolved to a uniform beam size of 
2'. Then the mean rms noise levels of the main beam 
brightness temperature were 0.18, 0.20, and 0.30K for 
the i^co (J=l-0), i^co (J=l-0), and ^^CO (J=2-l) 
lines, respectively. 

^^CO (J=l-0) emission data from the Boston 
University-Five College Radio Astronomy Galactic Ring 
Survey (BU-FCRAO GRS; Jackson et al. 2006) were also 
used. They offer a 46" angular resolution on a 22" grid 
and a 0.21 kms~^velocity resolution. 

2.2. Data of X-Ray and Other Wavebands 

The ROSAT X-ray image of the SNR Kes 78 region 
was gained from the ROSAT All Sky Survey broadband 

5http://www.iram.fr/IRAMFR/GILDAS 



data. The XMM-Newton data of partial region of Kes 78 
were obtained from an observation toward the Galac- 
tic ridge centered at (18''51'"44^67, 00°08'58".5, J2000), 
which was carried out on 2003 November 21-22 (Ob- 
sID: 0017740501, PI: F. Azita Vahnia) in the full frame 
mode and with the medium filter. The EPIC data were 
reprocessed with the Science Analysis System software 
(SAS, version 9.0.0). After removing time intervals with 
heavy proton flarings, net exposures of 29 ks and 27 ks re- 
mained for EPIC-MOS and EPIC-PN, respectively, and 
were used for analysis. The XMM-Newton field of view 
(FOV) covers only the northeastern boundary region of 
Kes 78, which is the brightest region in the 1.4GHz radio 
continuum emission. 

The Spitzer 24 fim mid-IR observation used here was 
carried out as 24 Micron Survey of the Inner Galac- 
tic Disk Program (PID: 20597; PI: S. Carey) with the 
Multiband Imaging Photometer. The 24 /im Post Basic 
Calibrated Data were obtained directly from the Spitzer 
archive. The 1.4 GHz radio continuum emission and the 
Hl-line data were obtained from the archival VLA Galac- 
tic Plane Survey (VGPS; Stil et al. 2006). 

3. RESULTS 
3.1. Molecular Environment 
3.1.1. The OH Maser Point and MCs at Vlsr ~ 81kms~^ 
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Fig. 1.— i^co (J=l-0) (solid line),i2cO (J=2-l) (dotted line), 
and ^^CO ( J=l— 0)(thick solid line) spectra at the maser point. 
The upper x-axis refers to the rest frequency of the ^■^CO (J=l-0) 
line. AH the spectra are not convolved to a uniform beam size. 



We first focus on the molecular gas indicated 
by the 86 kms^H720 MHz OH maser located at 
(18''51'°48".04,-00°10'35"). The CO spectra at the 
maser point are given in Figure [1] There are four promi- 
nent peaks at Ksr ~ 10, 37, 81, and 102kms-i. The 
LSR velocity of the maser, 86 kms~^, is in the red (right) 
wings of the Slkms^^ ^^CO lines, in which ^■^CO emis- 
sion is insignificant. A few secondary ^^CO components 
ranging from 52kms^^ to 96kms~^, with few ^"^CO 



3 




80 
V (km/s) 

Fig. 2. — Grid of CO spectra restricted to the velocity range 50kms~^ to 105kms~^ in the 1720 MHz OH maser region with 1' spacing 
(labeled as the box "M" in the bottom right panel of Figure (3)1. Note that the ^^CO (J=2— 1) lines have been multiplied by a factor of 
1.6 so as to scale the ~ 81kms~^ ^■^CO (J=2— 1) peak of the maser point ("Me") with the ^^CO (J=l-0) peak. All the spectra are not 
convolved to a uniform beam size. 



counterparts, seem to be linked with the prominent com- 
ponent peaked at 81kms~^. As ^^CO emission, usuahy 
optically thin, arises in the quiescent dense gas along the 
whole LOS, the lack of significant ^'^CO features at the 
corresponding LSR velocities indicates that these ^^CO 
secondaries are likely to represent the disturbed gas devi- 
ating from the systemic velocity ~ 81 kms""'^, as compli- 
cated broadened wings. Similar CO line profiles are also 
present at points close to the OH maser (point "Me"), as 
seen in the grid (Figure [2|). By scaling the ~ Slkms"^ 
i^CO (J=2-l) peaks to match the ^^CO (J=l-0) peaks, 
with the maser ("Me") as the fiducial point, relative en- 
hancement of the ^^CO (J=2-l) is seen at some points 
(e.g., "Mb", "Mc", and "Mf at ~ 90kms-i and "Ma", 
"Mc", and "Me" at ~ 67km s^^), which may be indica- 
tive of relatively warm shocked gas. The shock interac- 
tion is also favored by the elevated ^^CO J=2-l/J=l-0 
ratio at 64-68 kms~^ in the maser vicinity (see below). 

By inspecting the intensity maps in the entire veloc- 
ity range, however, we only found two velocity intervals 
(around 67kms~^ and 81kms~^) in which the ^^CO 
emission demonstrates more or less morphological corre- 
spondence with the SNR. As seen in Figure [Sj at round 
81 kms~^, the CO emission exhibits a cavity structure, 
which is consistent with the SNR as delineated by the ra- 
dio contours. The eastern bright radio shell of the SNR 
overlaps with most of the region of strong ^^CO emis- 
sion, which covers the maser point. The radio emission 
fades from east to west, and a ^^CO arc composed of a 
few clumps is located along the faint radio boundary in 
the west. In the SNR boundary, four prominent ROSAT 
0.1-2.4 keV X-ray patches are seen, two of which are lo- 
cated along the western radio shell and one of which is 
coincident with one of the western ^^CO clumps ("clump 
1", see Section 13X21). 

At KsR ~ 67kms-^ (Figure O, a string-like ^^CO 
feature extends from the northeast to the south in the 
FOV, also covering the maser point, and most of it seems 



to follow the western radio shell. However, two bright 
ends of the string are located outside the remnant, and 
little ^^CO emission corresponds to the radio brightness 
peak on the northern shell. 

The integrated CO-line intensity ratio between the dif- 
ferent transitions is used as kinematic evidence of inter- 
action between SNR shocks and MCs (e.g.. Seta et al. 
1998; Jiang et al. 2010). Elevated ^^CO J=2-l/J=l-0 
line ratios ~ 1 are also present on the boundary in the 
two velocity intervals of concern. At 76-79 kms-\ the 
^^CO J=2-l/J=l-0 ratios are apparently elevated to 
> 1 (compared with ^ 0.4-0.6 in average) at several po- 
sitions along the radio boundary. Apart from the north- 
western and southern boundaries, where an incomplete 
grid of points spaced every 2' was observed in PMOD for 

CO ( J=l-0) and thus the ratios there should be taken 
with caution, the ratio elevation in the northeast and 
the west is noteworthy. The northeastern ratio elevation 
(w 1) at around (18h51'"37", -00°01'22", J2000) is lo- 
cated in a bright patch of radio emission, and the west- 
ern elevation (> 0.9) around (18''50'"16", -00°13'00", 
J2000) is coincident with a molecular clump ("clump 
1", see Section r3.1.2p . In some interacting SNRs, the 
line ratios are observed to exceed unity (see Dubner et 
al. 2004), as found in SNR 3C396, however, due to some 
reasons such as beam dilution of the small filling factor 
of the ^^CO (J=2-l) emitting region compared to that 
for ^^CO (J=l-0), the observed relatively elevated ra- 
tios > 0.9 can also be a good probe of the SNR-MC 
interaction (Su et al. 2011). At 64-68 kms^"'^, the line 
ratio is enhanced to > 1 in a region around (18'^51™55*', 
-00°10'00", J2000) adjacent to the OH maser point. By 
comparison, the line ratios are low in the two bright 
ends, which are outside the SNR boundary, of the east- 
ern 67kms~^ ^^CO string. The high line ratio near the 
maser point is another signature of disturbance of the 
molecular gas there by the SNR shocks and the gas emis- 
sion at around 67kms~^ is in the blue wings of the main- 
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TABLE 1 

Some Parameters of the Molecular Gas Associated with 
Kes 78 



Regions 


Ar(H2)(10 


2icm-2) 


M{Mq) 




X 


L'i'E 


X L'i'E 


Maser point 


8.8 


8.5 




Clump 1 


4.3 


3.0 


1.9 X 10^ d§ 1.3 X 10^ dl 



Note. — Here two methods (indicated as "X" and "LTE", re- 
spectively) are used to estimate t he par ameters of the molecular gas; 
see the last paragraph in Section l3.1.1l 

body clouds at 81kms^^. 



Figure H shows a large FOV of the ^^CO (J=l-0) en- 
vironment at 80-84 km s"^ around SNR Kes 78. The 
SNR is seen to be located in the side of some MCs. A 
long molecular strip in the northeast seems to connect 
the radio shell, although the LOS projection effect can- 
not be excluded, and the radio brightness peak along the 
northwestern SNR boundary is located at the connec- 
tion point. The SNR also seems to be in contact with a 
molecular patch in the east and the OH maser is located 
at the interface. The eastern molecular patch seems to 
spatially correspond to some extent to the extended VHE 
7-ray source HESS J1852-000 (considering the relatively 
large point-spread function of the H.E.S.S. observation) .31 
The ^'^CO emission is much fainter in the west of the rem- 
nant than in the east. This, together with the faint radio 
emission in the west of the remnant, is consistent with 
a scenario that the remnant blows out to a low-density 
region. We also inspected the ^^CO (J—l-O) emission 
at around 67kms~^, but no distinct molecular structure 
was seen. 

The above evidence from CO emission, such as the 
morphological correspondence of the CO emission with 
the SNR, the enhanced ^^CO J=2-l/J=l-0 ratio, and 
the ^^CO secondary peaks around 81kms~^ as compli- 
cated broadened features, suggests that SNR Kes 78 is 
associated with the molecular gas at a systemic velocity 
of 81 kms~^. There is a discrepancy of ^ 5kms^^ be- 
tween the LSR velocity of the OH maser and the systemic 
velocity of the associated MCs. A similar discrepancy is 
also seen in SNR W28, in which the velocity offset be- 
tween the OH masers and the main-body MC is as high 
as 8kms^^ (Claussen et al. 1997). 

The column density of the molecular gas, A^(H2), at 
the OH maser point (within the beam size 56" of the 
PMOD observation) is estimated with two methods and 
presented in Table [T] In the first method, A^(H2) is cal- 
culated by using the ^^CO (J=l-0) velocity-integrated 
brightness temperature together with the mean CO-to- 
H2 mass conversion factor A^(H2)/PF(^^C0) w 1.8 x 
10^° cm-^K-^km^^s (Dame et al. 2001). In the sec- 
ond method, on the assumption of local thermodynam- 
ical equilibrium (LTE) with an excitation temperature 



15 K and ^^CO ( J=l-0) being optically thick, the "CO 
column density is converted to A^(H2) using the relation 
A^(H2) « 7 X 10^ A^("CO) (Frerking et al. 1982). The 
results obtained from the two methods are similar. 

3.1.2. The Western Molecular Arc 

A close-up tri-color image of the western molecular arc 
is presented in Figure [5ja). As previously mentioned, an 
X-ray patch is coincident with "clump 1" along the arc. 
The X-ray patch seems to extend northwards to another 
small clump ( "clump 2" ). (The two clumps are labeled in 
Figure [H) A 24 /im mid-IR patch also seems coincident 
with "clump 1". Another large, bright mid-IR patch in 
the north part of the image appears to cover the north- 
ernmost CO clump and, projectively, an ultracompact 
HII region at the velocity range of 7-26 km s~^ (also see 
Section I4.4.2|) . Interestingly, some points with elevated 
12C0 J=2-l/J=l-0 ratios (> 0.9) at 76-79 kms"! are 
seen at and around "clump 1" (see Figure[3l bottom right 
panel). Such ratios significantly higher than the average 
(0.4-0.6), as mentioned in Section [3.1.11 are also a sig- 
nature of warm shocked molecular gas. 

For the cores of two of the molecular clumps, two 
grids of CO spectra are presented in Figures |6] and [3 in 
which broad ^^CO line profiles at around 81kms~^ can 
be found. For "clump 1", blueshifted (left) broadenings 
of the ^^COline profiles are seen at points "la" (extend- 
ing to at least 73kms-i), "lb", "le", and "Ih". Bumps 
in the red wings of ^^CO lines are seen at "lb" , "Ic" , and 
"If without corresponding ^'^CO emission. The broad 
red wing of ^^CO (J=l-0) at "le" is not decided be- 
cause the ^'^CO emission is not ignorable {> 3a) in the 
86-90 km s~^ interval. For "clump 2", blueward broad- 
ening of the ^^CO line wings is seen at "2a", "2c", and 
"2f ' , while a redward broadening is clearly seen at "2d" . 
These ^^COline broadenings in red or/and blue side(s) 
are as large as about 10kms~^, which is strong evidence 
that the two clumps are shocked. 

3.1.3. Other Molecular Components along the Line of Sight 

We also find some noteworthy molecular structures at 
different LSR velocities in the Kes 78 field. Figure [8] 
shows the Ksr ^ 100 kms^^ ^^CO (J=l-0) map com- 
pared with the Spitzer 24 fim mid-IR emission image. 
The CO emission at this LSR velocity (the velocity at 
the tangent point) is strong in the northeast region of 
the map (projectively across the northeastern border of 
the SNR) and, interestingly, is similar to the distribution 
of the bright 24 fim diffuse emission (also note a hollow, 
in the both wavelengths, outsides the remnant). In both 
of the ^■^CO and 24 /im emission maps, there is a north- 
south oriented branch stretching from the northeastern 
bright portion and well following the radio continuum 
"ridge" that is usually thought to be a structure of the 
remnant. At this LSR velocity, however, no molecular 
features other than the "branch" correspond to the other 
structures of the SNR (e.g., especially, the eastern radio 
shell). As will be discussed in Section l473l the faint radio 
ridge is likely to be irrelevant to the SNR. 

Serendipitously, in the LSR velocity interval 86- 
97kms~^, we discern a bubble/ring-like ^^CO (J=l-0) 
structure of a radius ^ 15' (see Figure [9l) centered at 
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Fig. 3. — Upper and middle rows: intensity maps of ^^CO (J=l-0) and ^^CO (J=2-l) at the Vlsr ~ 67kms~^ and Vlsr 81 kms"^. 
Thick contours (red) represent the ROSAT 0.1-2 keV X-ray emission with levels of 0.22, 0.26, and 0.30 counts pixel"'^ (after being smoothed 
with a 3' Guassian kernel). Bottom row: ^^COJ=2-1/J=1-0 ratio maps in 64-68 kms-l and 76-79 kms"!, after convolving the PMOD 
l^CO (J=l-0) to the same beam size as the KOSMA ^^CO (J=2-l) data. The thick contours (green) represent the ^^CO (J=2-l) 
emission at 80-84 kms~l interval with levels 6.5, 10.1, and 13.7 K kms"^ Here both ^^CO{J=2-l) and ^^CO(J=1-0) data used achieve 
a signal-to-noise ratio > 3. Thin contours (black) of VGPS 1.4GHz continuum emission are overlaid in six panels with levels of 16, 20, 24, 
28 and 32K. The plus sign in each panel denotes the OH maser point. The boxes, labeled with "M", "1", and "2", define the OH maser 
region, the cores of "clump 1", and "clump 2", the spectra of which are shown in Figures [2l l6l and [71 



about (18''51'"30", -00°09'00", J2000). Although it ap- 
pears in almost the same direction as SNR Kes 78, with 
am angular size similar to that of the SNR, they are 
not completely coincident with one another and there is 
no corresponding feature between them. The systemic 
LSR velocity 91kms~^ indicates that the bubble/ring is 



located at the near distance ~ 5.5 kpc or the far dis- 
tance ^ 7.9 kpc, in the background of the SNR (see 
Section 13.21 for the method for deriving the kinematic 
distance). Using the second method for estimating the 
molecular column density described in Section FS.l.ip . the 
mass of the bubble-like molecular gas is estimated to be 
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Fig. 4. — Molecular environment of Kos 78 in the 80-84 km s ^ 
interval (the BU-FCRAO GRS ^^CO (J=l-0) observation), over- 
laid with 1.4 GHz radio contours at levels 16, 19.75, 23.5, 27.25, and 
31 K. The plus sign denotes the OH masor point. The two boxes 
indicate the spatial ranges of the two maps given in Figure [5] 

~ 3.1 X 10^ Mq and ~ 6.4 x 10^ Mq for the near and 
far distances, respectively. The n ature of the structure 
will be discussed in Section 14.4.31 

3.2. The kinematic distance to SNR Kes 78 

Due to the association of Kes 78 with the MCs at 
the systemic LSR velocity ~ 81kms~^ (Section I3.1.ip . 
there are two candidate kinematic distances to the SNR, 
4.8 kpc (near side) and 8.6 kpc (far side). Here we 
have used the Clemens' (1985) rotation curve of the 
Milky Way together with i?o = 8.0 kpc (Reid 1993) and 
Vb = 220 km s-^ The LSR velocity of the OH maser, 
86kms~^, corresponds to similar distances 5.2 kpc and 
8.3 kpc. 

A modified HI absorption method (Tian et al. 2007) 
is next used to discriminate between the near and far 
distances. The HI spectra are extracted from two radio- 
bright regions as labeled in Figure |8] (regions "I" and 
"H" ) and plotted as shown in Figure [101 A few dis- 
tinct (> 3(t) absorption features appear below Vlsr 
100 kms^^, such as at lOkms"^ and in the 80-90 kms^^ 
interval. However, there is no absorption between 
95kms~^and llOkms^^, which corresponds to the tan- 
gent point (at a distance around 6.7 kpc). Therefore, 
SNR Kes 78 is in front of the tangent point, namely at 
a distance smaller than 6.7 kpc; hence, the remnant is 
at the near distance 4.8 kpc. Hereafter we parameterize 
the distance as d = 5^5 kpc (also considering the LSR 
velocity discrepancy of the OH maser from the systemic 
velocity) . 



3.3. XMM-Newton X-Rays from the Northeastern 
Boundary 

3.3.1. X-Ray Image 

From the available XMM-Newton EPIC-MOS X-ray 
data, we produce a tri-color X-ray image of the north- 
eastern boundary of Kes 78 (as shown in Figure [5l panel 
(b); 0.5-1.0 keV coded in red, 1.0-2.0 keV in green, and 
2.0-7.0 keV in blue), which is spatially correspondent 
to the radio brightest region. In the image production, 
the background subtraction and exposure correction have 
been applied, and adaptive smoothing has been made 
with a signal-to-noise ratio (S/N) ratio > 3. 

In the image, the X-ray-emitting region is well bor- 
dered by the radio contours, which indicate a sharp radio 
brightness drop at the blast shock. This part of X-ray 
emission was not detected in the ROSAT X-ray survey. 
The eastern portion of the shown region (dominated by 
the orange color) is apparently softer than the western 
portion (dominated by the green color), which is coinci- 
dent with the radio brightness peak. 

3.3.2. Spectral Analysis 

Before the spectrum extraction, point-like X-ray 
sources in the FOV were detected with the wavelet 
method and removed from the events file. We then ex- 
tracted the on-SNR X-ray spectra from two regions (la- 
beled with "A" and "B" , as shown in Figure [Til the 
raw MOS image), which cover the hard (green) and soft 
(orange) portions (mentioned above), respectively. Be- 
cause of the small amount of available counts, the MOSl 
and M0S2 spectra were merged to increase the statis- 
tical quality. From the EPIC-PN observation, only the 
spectrum of "region A" was obtained. Thus for region A, 
both the MOS and PN spectra were used simultaneously, 
while for region B, only the MOS spectrum was used. 
As the X-ray source region is located near the edge of 
the FOV and is not very bright relative to the instru- 
mental background, a double background subtraction 
method was used. The off-SNR source-free spectra were 
extracted from both the MOS and PN data (the extrac- 
tion regions, as also shown in Figure lll[ do not include 
the significant X-ray emission from the Galactic ridge 
in the center of the FOV). The respective instrumen- 
tal background contributions from the same on- and off- 
SNR regions were obtained from the Filter Wheel Closed 
data of MOS and PN and were subtracted from the on- 
and off-SNR spectra. The individual on- and ofF-SNR 
spectra were adaptively binned to achieve a background- 
subtracted S/N of 3. 



The XSPEC spectral fitting package (ver.11.3) was 
used. Each on-SNR spectrum was jointly fitted together 
with the off-SNR spectrum (Figure [12). The on-SNR 
diffuse background was determined by scaling the off- 
SNR emission according to the region sizes (Figure [I2\i 
and is phenomenologically well described by an absorbed 
nei + blackbody model. For the foreground absorption, 
the cross sections from Morrison & McCammon (1983) 

^ http: / /heasarc.gsfc. nasa. gov /docs/ software /Iheasoft / xanadu / 
xspec/xspecll/index.html 
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Fig. 5. — (a) A close-up tri-color image of the western molecular arc (the region of the image is indicated by the solid box in Figure|4]l. 
Red: Spitzer 24 fira emission; green: intensity map of ^^CO (J=2-l) integrated from 80 kms~^ to 84 km s"'^ (overlaid with contours at 
levels of 6.5, 8.3, 10.1, 11.9, and 13.7 K kms"!); and blue: ROSAT 0.1-2.0 keV X-ray map (smoothed to 3'). (b) Tri-color XMM MOS 
X-ray image of the northeastern part of Kes 78 (the region of the image is indicated by the dashed box in Figure |4]l. The X-ray intensities 
in the 0.5—1.0, 1.0—2.0, and 2.0-8.0 keV bands are coded in red, green, and blue, respectively. The 1.4 GHz radio contours are at the same 
levels as those in Figure |3] 
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with 1' spacing (as labeled in the 



were used, and solar abundances were assumed. The 
spectra of both regions A and B show distinct line fea- 
tures of Mg Hea (~ 1.35 keV) and Si Hea (~ 1.85 keV), 
indicating the thermal origin of the emission. The net 
SNR X-ray emission can be well fitted with an absorbed, 
non-equilibrium ionization thermal plasma model ( vnei) 
with solar abundances. The fitting results are summa- 
rized in Table [H The physical properties between the 
two regions are similar, except that the intervening hy- 
drogen column density for region A is higher than that 
for region B. The relatively high A'h for region A is con- 
sistent with a contact of the SNR at the northeastern 
border (coincident with region A) with the dense molec- 



ular strip (the shock strip interaction can naturally ex- 
plain the radio brightness peak there; see Section 13.1.11 
and Figure |4]) and explains the fact that X-rays from re- 
gion A are harder than those from region B (see Section 
I3.3.1|) . The gas temperatures are found to be 1.0-1.9 keV 
for region A and 0.7-2.9 keV for region B. The ionization 
timescales Ugti of the X-ray-emitting gas of the two re- 
gions are both on the order of 10^° cm~''s, much smaller 
than 10^^ cm~'^s, which indicates that the hot gas has 
not yet reached ionization equilibrium. Estimates of the 
hydrogen number densities, tih, derived from the volume 
emission measures {fuennV, where / is the filling fac- 
tor of the hot gas and rig ^ 1.2nH is assumed) are also 
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TABLE 2 

Spectral Fitting Results for the Northern Shell of Kes 78 with 
90% Confidence Ranges and Estimates of the Gas Density 



Regions 


A 


B 


Net count rate (10~^ counts s^^) 

Xl (dof) 

A^H (10^2 cm-2) 

kT^ (keV) 

n^U (10l°cm-3s) 

frieUiiV/dl (1056cm-3) 

Flux(10-12 ergcm-2 s'^Y 


3.61 ±0.11 (MOS) 
8.52 ±0.27 (PN) 
1.04 (318) 

i-'J^-o.os 

1 ^1+0.34 

-I 09+0.89 
-'-■°^-0.32 
1 -,,+0.31 
-'-■-'-^-0.25 
o 41+0.94 

ii+o-oi 

'-'■^^-0.01 


2.79 ±0.10 (MOS) 

1.17 (137) 
n 70+012 

'J- "^-0.09 

1 sn+1'5^ 

l-'^'J-0.65 

0.58l°i| 

1 or, + 2.19 

l-°^-0.92 
'^•^'-'-0.02 



Note. 



The net count rates listed in the table are for the on-SNR 



spectra. The rates of the off-SNR spectra are 6.17 ± 0.18 X 10 counts s 
and 1.26 ±0.04 X 10"^ counts s"! for MOS and PN, respectively. 

The unabsorbed fluxes are in the 0.6-5.0 keV band. 
^ In the estimate of the densities, we assume oblate spheroids for elliptical 
regions A (with half-axes 3. '36 X 3.'36 X 1.'84) and B (2. '74 X 2.'74 X 1.'75). 
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Fig. 7. — Grid of CO spectra restricted to the velocity range 
65kms~i to 95kms~i in "clump 2" with 1' spacing (as labeled 
in the bottom right panel of Figure |3j . All the spectra are not 
convolved to a uniform beam size. 

given in Table [2j In the derivation of the densities we 
have assumed that the three-dimensional shapes of the 
elliptical regions (A and B) are ablate spheroids. Devia- 
tions from these assumptions and the non-uniformity of 
the X-ray-emitting plasma are consolidated into factor / 
of the individual regions. 

4. DISCUSSION 

4.1. Global Evolution of the SNR and the Progenitor 

As revealed above (Section [XTHl), SNR Kes 78 evolves 
in an interstellar environment of molecular gas with sys- 
temic velocity ~ 81 kms~^ and expands into a cloudy re- 
gion in the east. Interestingly, the X-ray-emitting gas in 
the northeastern edge is characterized by a single ther- 
mal temperature (around 1.5 keV) and a low number 

density {n^ - Q.lf-^/'^d^^''^ cm This indicates that 
the blast wave in the northeast of the SNR propagates 
in the intercloud/interclump medium (ICM). 



The velocity of the blast wave propagating in the 
ICM can be estimated as Vg — [l&kTx/ {'ip.rau)Y^'^ 
1-1-0 2 ^ lO'^kms"^, where the fitted hot gas tempera- 
ture for region A (for which the error bars of the spec- 
trally fitted parameters are smaller than those for region 
B), kTx = 1.5llg;gg keV, is adopted as the postshock gas 
temperature and the mean atomic weight p. equals 0.61 
for fully ionized plasma. We take the curvature radius of 
the eastern radio shell ^ 12' as the SNR radius, namely, 
rs ~ 17c?5 pc, assuming that the dynamical age of the 
remnant, t ~ 2rs/{5vs) ~ 6.1^q ^ kyr, is derived with 
the Sedov (1959) evolution law. The explosion energy 
is £; = (25/4C)(1.4?iomH)r3i;2 _ 4.7 x 10^^° f-'^/^dl^^ erg, 
where ^ — 2.026 and uq — 71h/4 is the prcshock inter- 
cloud hydrogen density. The ionization age of the X-ray 
emitting gas in the northeastern edge is inferred from 
the ionization timescales (see Table for region A, 



t,. 



4.4^0 9 kyr old, which is slightly, but reasonably, 
smaller than the dynamical age. Both estimates suggest 
a relatively young age of Kes 78. 

The absence of [OIII] leads to an inference that the 
velocity of the shock wave responsible for the optical fil- 
amentary and diffuse structures is < lOOkms"^ (Boumis 
et al. 2009). This shock cannot be the SNR blast wave, 
otherwise the SNR should have been in the radiative 
phase with a considerably large age (Chevalier 1974): 
0.31rs/vs ^ 52 kyr, inconsistent with the X-ray proper- 
ties obtained above. As Boumis et al. (2009) suggest, 
however, it is in the interstellar clouds that the slow 
shocks are propagating. 

The bright radio emission in the east of the SNR is 
in agreement with the aggregation of the ^ 81kms~^ 
molecular gas there and the inference that the clouds are 
clumpy. The radio continuum emission may arise from 
the blast shock that propagates in the ICM (Blandford 
& Cowie 1982). Contrarily, in the west, the ~ 81kms~"'^ 
CO emission is weak and the radio continuum fades out. 

As also mentioned in Section 13.1.11 in the west, the 
SNR expands into a relatively low-density region; but 
there the SNR seems to collide with a clumpy molecular 
arc. The entire SNR appears to be in a cavity surrounded 
by molecular matter. The molecular number density of 
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Fig. 8. — Spitzer 24 /xm map (left panel) and BU-FCRAO GRS ^^CO (J=l-0) antenna temperature map integrated from 99 kms ^to 
101 kms~^(right panel ). Both of the maps are overlaid with 1.4 GHz radio emission contours with the same levels as those in Figure [3l 
The pl us si gn in each panel denotes the OH maser point. The two boxes lab eled with "I" and "H" are regions from which the HI spectra 
(Figure (ToJ are extracted for discriminating two candidate distances (Section l3.2|l . 




Fig. 9. — Image of the '^^CO (J=l-0) bubble-Uke structure at 
86-97 km s-i discerned from the BU-FCRAO GRS data, overlaid 
with 1.4 GHz radio emission contours with the same levels as those 
in Figure |4] The plus sign denotes the OH maser point. 



"clump 1" along the arc is ' 



270d^^ cm' 



7188^5 



cm 




with an angular radius 2' adopted. The broad CO line 
wings and enhanced ^^CO J=2-l/J=l-0 ratios along 
the western arc impliy a disturbance suffered by the 
molecular gas. Similar molecular arcs are seen in some 
other interacting SNRs, such as Kes 69 (Zhou et al. 2009) 
and Kes 75 (Su et al. 2009), in which the arcs are sug- 
gested to represent the debris of the clumpy shells of 
the progenitors' wind bubbles. If the molecular cavity in 
Kes 78 was created by the progenitor, using the linear re- 
lation between the main-sequence wind bubble size and 
the progenitor's mass (Y. Chen, et al.. 2011, in prepa- 
ration Q), i?b = I.2IM/M0 - 8.98 pc, the cavity radius 

^http: / /astronomy.nju. edu.cn / ~ygchen / papers /Rb-M /Rb- 
M.pdf 



LSR (km/s) 



Fig. 10. — VGPS HI spectra extracted from the regions "I" and 
"II." The long-dashed lines represent the 3(t level. 



17pc would imply a mass around 22Mq for the progeni- 
tor star of SNR Kes 78, which seems to be of 09 type. 

4.2. Pressures in the multi-phase gases 

We have arrived at such a scenario in Section |4?11 that 
SNR Kes 78 evolves in a cloudy interstellar medium, 
the blast wave in the ICM giving rise to the X-rays 
and the cloud shocks in the dense clouds to the optical 
emission. The thermal pressure of the X-ray-emitting 

gas is 2.3nHT^ ~ 6.0 x lO^f-^^^d^^^^cm-^K and ~ 
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Fig. 11. — XMM-Newton raw image of Kes 78. Two solid el- 
lipses (labeled "A" and "B") and three dotted boxes are defined, 
respectively, for source and background spectrum extraction. The 
contours of 1.4GHz radio emission outline the northeastern area of 
Kes 78, which have the same levels as those in Figure |4] 

4.0 X lO^/^^/^rfjT^^^ cm~^K for regions A and B, respec- 
tively. According to the density-sensitive line ratio of 
[SIX] AA6716/6731, B09 determined the electron densities 
in the shocked clouds to lie below 200 cm^^ and suggest a 
temperature of the electrons ~ 10"* K. The upper limit of 
the thermal pressure in the optically emitting structures 
is estimated to be 4 x 10^ cm~^K, which is similar to 
the pressure derived for the X-ray-emitting gas. Thus, it 
seems that the shocked ICM may be in pressure balance 
with the shocked clouds. 

However, there is significant pressure variation in the 
multi-phase cloudy gases in the east of the remnant. The 
1720 MHz OH masers are believed to arise from very 
dense (~ 10^ cm~^) regions of MCs that are shocked 
by the C-type shock waves (Lockett et al. 1999). By 
means of Zeeman splitting, the LOS component of the 
magnetic field in the Kes 78 OH maser region, -Bios, was 
measured to be 1.5±0.3mG (Koralesky et al. 1998). For 
a randomly oriented field, the post shock field strength 
Bp, = 2Bios (Frail & Mitchell 1998; Crutcher 1999). 
Thus, the magnetic pressure in the OH maser region of 
Kes 78 is BlJiSnk) w 2.6 x lO^cm'^K, which is over 
two orders of magnitudes higher than the thermal pres- 
sures of the X-ray emitting intercloud plasma and the 
optically emitting filamentary and diffuse structures. No- 
tably, similar pressure variation from the X-ray emitting 
gas to the OH maser regions is seen in SNR W28 (Rho 
& Borkowski 2002). Such a pressure contrast has been 
suggested to be caused by the interaction of the radiative 
shell with the dense molecular clumps (Chevalier 1999). 

4.3. The Dubious Ridge 

The radio "ridge" was sometimes suspected to be the 
western part of an elongated SNR shell (e.g., Caswell et 
al. 1975). We have seen in Section 13.1.31 that the radio 
ridge appears to be coincident with the northsouth ori- 
ented 24 fiia IR branch, which extends from the bright 
emission in the northeastern region outside the remnant's 
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Fig. 12. — XMM-Newton EPIC spectra of regions A and B 
defined in Figure 1111 The on-SNR spectra (black solid lines for 
EPIC-MOS and green solid line for EPIC-PN) are fitted as a sum 
of the diffuse background emission and the net SNR emission. The 
on-SNR diffuse backgrounds (red dotted lines for MOS and blue 
dotted line for PN) are scaled from the ofT-SNR spectra (red solid 
lines for MOS and blue solid line for PN), which are mimicked by 
an absorbed nei + blackbody model. The net SNR spectra (black 
dashed lines for MOS and green dashed line for PN) are fitted well 
with an absorbed nei model. 

radio boundary, the 24 /im emission in and around the 
eastern region of the SNR has a strikingly similar dis- 
tribution to the KsR ^ lOOkms-i "CO emission. This 
leads us to suspect that the radio ridge is not a part of the 
SNR that is associated with the ~ 81kms"^ MCs, but 
rather is associated with the ~ 100 kms^^ MCs, which 
is essentially located at the tangent point of the inner 
Scutum-Centaurus spiral arm of the Galaxy. On the 
other hand, optical observation towards Kes 78 finds no 
significant [SII] from the radio ridge, while in the eastern, 
northern, and southern regions of Kes 78, the [SII] emis- 
sion is evident, with the [SII]/Ha ratio > 1.3 (Boumis 
et al. 2009). As the [SII] strength in the shock- ionized 
gas is comparable to the Ha strength and is much fainter 
in the photoionized gas (HII regions or planetary nebula; 
Dennefeld & Kunth 1981), the radio emission is probably 
of thermal origin and possibly arises from star formation 
region(s). Furthur study estimating the radio spectral 
index for the radio "ridge" will test its nature. 

4.4. Other Sources in the Kes 78 Field 

AAA. PSR1850-0006 

Since the progenitor of Kes 78 was probably a massive 
star (see Section |4T|) . a neutron star or pulsar may have 
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been left after its explosion. PSR1850-0006 is projected 
in the northwestern portion of the SNR (Keith et al. 
2009). However, with a rotational period of 2.2 s and a 
period derivative of 4.3 x 10~^^ ss~^, the pulsar has a 
characteristic age of about 8 Myr, which is much older 
than the age of SNR Kes 78 (Section l4.ip . disfavoring an 
association between PSR1850-0006 and the SNR. 

4.4.2. HII Regions 

The uhracompact HII region located at (18''51'°25^ 
00°04'07", J2000) is northeast of the SNR, as seen in 
Figure |4l It is coincident with the massive star forma- 
tion region IRAS 18488+0000 (Bronfman et al. 1996) and 
may be associated with SNR Kes 78. In fact, a 22 GHz 
H2O maser at Vlsr 80kms^^, which is very similar 
to the systemic LSR velocity of the MCs with which the 
SNR is in contact, is suggested to be associated with the 
star formation region and is at the same distance (5 kpc) 
as we determine here for the SNR (Beuther et al. 2002). 
The HII region or IRAS 18488+0000 is located at the 
northern molecular strip with which the SNR interacts 
and only 2' north of the radio brightness peak of the 
SNR. 

Another ultracompact HII region at (18''50™31^ 
-00°01'55", J2000) (Garay et al. 1993) appears to be 
projectionally located on the northwestern boundary of 
the SNR and connects the western molecular arc in the 
northwest, as seen in Figure |H In this HII region, six 
22 GHz H2O masers were detected at the LSR velocity 
range ~ 7-26 kms"^ (Hofner & Churchwell 1996). There 
may not be associated with the SNR because they are at 
different LSR velocities. 

4.4.3. The Molecular Structure at Vlsr ~ 91kms~^ 

As concluded in Section I3.1.3| the bubble/ring-like 
structure discovered at Ksr = 86-97 kms~^ is at a 
furthur distance (5.5/7.9 kpc) than that to SNR Kes 78 
(4.8 kpc). If it is an interstellar bubble, then the mean 
molecular number density of the undisturbed gas before 
the bubble was created is ~ 77/54 cm~^ for the near/far 
distance. Adopting an expansion velocity Vb ^ 7kms^^ 
(FWHM of the ^^CO line at 91 kms-^), the kinematic 
energy is ~ 1.5/3.1 x 10^° erg. If the bubble is blown by 
a wind , it has an age 3rb/{5vb) ^ 2.0/2.9 x 10^ yr and 
the mechanical luminosity of the wind is ~ 1.2/1.8 x 
]^Q37gj.gg-i according to the canonical evolutionary law 
of a wind bubble given by Weaver et al. (1977). This 
ring/bubble-like structure may not be a distinct SNR. 
Actually, for an SNR in the radiative stage (Chevalier 
1974), its age is 0.31r,,/?;f, ~ 1.0/1.5 x lO^yr, but the 
supernova explosion energy would be 4.6/9.6 x 10^^ erg, 
which is much higher than the fiducial value 1 x 10^^ erg. 

5. SUMMARY 

We have investigated the molecular environment of 
the Galactic SNR Kes 78 using the CO observations in 
PMOD, KOSMA, and FCRAO and have performed an 
XMM-Newton X-ray spectroscopic study for the north- 
eastern edge of the remnant. The main results and con- 
clusions are summarized as follows: 



1. SNR Kes 78 is found to be associated with the MCs 
at a systemic LSR velocity of 81 kms~^. At around 
this velocity, the SNR is revealed by the ^'^CO ob- 
servation to be in the side of some dense MCs and 
is consistent in extent with a cavity of ^^CO gas. 
The SNR expands into a cloudy dense region in 
the east, where the OH maser emission arises and 
the bright radio shell overlaps with most of the re- 
gion of strong ^^CO emission. Broadened ^^CO 
line profiles discerned in the eastern maser region 
and the western clumpy molecular arc, the elevated 
i^CO J=2-l/ J=l-0 ratios along the SNR bound- 
ary may be signatures of shock perturbation in the 
molecular gas. 

2. The association of SNR Kes 78 with the ^ 
81kms^^ MC, together with the HI absorption 
along the LOS, places the SNR at a kinematic dis- 
tance of 4.8 kpc. 

3. The X-rays arising from the northeastern radio 
shell are emitted by underionized hot 1.5 keV), 
low-density (^ 0.1 cm"'^) plasma with solar abun- 
dance, and the plasma may be of intercloud ori- 
gin. Especially, the X-rays from the radio bright- 
ness peak, where a long molecular strip appears 
to connect, suffers a heavier absorption than those 
from other parts along this XMM-observed section 
of the shell. The thermal pressures of the X-ray- 
emitting gas and the optically emitting structures 
along the shell are lower than the magnetic pres- 
sure in the OH maser region by up to two orders 
of magnitudes. The age of the remnant is inferred 
to be about 6 kyr. 

4. The size of the molecular cavity in Kes 78, which is 
assumed to be created by the main-sequence wind, 
implies the progenitor's mass to be around 22M0. 

5. The north-south oriented radio continuum ridge 
west of the bright shell is suggested to be irrele- 
vant with the SNR, but related to a ^ lOOkms"^ 
CO branch at the tangent point, which corresponds 
to IR bright, low-[SII]/Ha ratio star formation re- 
gion(s). 

6. Along the LOS toward SNR Kes 78, we discern an 
irrelevant bubble/ring- like ^^CO structure 15' in 
angular radius at the LSR velocity 91kms~^, but 
the possibility of it being a separated SNR is ruled 
out. 
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